Studies on the Effect of Rain, Fog, and Smog on 
Outdoor Optical Wireless Links 


A Thesis Submitted 

in Partial Fulfillment of the Requirements 
for the Degree of 

Master of Technology 


by 

MANISHA AGRAWAL 



Department of Electrical Engineering, 
Indian Institute of Technology Kanpur 

May, 2005 



n 


CERTIFICATE 


This is to certify that the work contained in the thesis entitled “Studies on the Effect of 
Rain, Fog, and Smog on Outdoor Optical Wireless Links”, by Manisha Agrawal 
(Roll No. Y3104051), has been carried out under my supervision and this work has not 
been submitted elsewhere for a degree. 


May, 2005 Dr. Joseph John, 

Professor, 

Department of Electrical Engineering, 
Indian Institute of Technology Kanpur, 
Kanpur -208 016. 



■f 2 JU!, 7m/£l' 

•TRcftW Sft^it^T^ 9>«rTfr «t»T^ 

«i».*.*...lji2ili£_I- 


aIS2 





Ill 


ABSTRACT 

An optical wireless system is an alternate wireless system for high speed data 
transfer to combat the highly congested RF spectrum. Outdoor optical wireless systems 
are becoming more popular and gaining market acceptance as a functional wireless tool, 
because of their advantages such as, unlimited and unregulated spectrum, low cost etc. 
One barrier which still exists in the wide acceptance of optical wireless communication 
(OWC) Systems is the effect of the atmospheric conditions viz. rain, fog, haze, smog etc. 

The main emphasis in this work is on the degradations induced by the atmosphere 
on the received optical signal. A comprehensive study of the existing outdoor optical 
links is done and their major features are compared. The major design challenges faced 
by an Outdoor wireless link are also discussed. The major atmospheric degradations 
considered are: attenuation due to fog, attenuation due to rain, and attenuation due to 
smog. Theoretical models for rain, fog and smog are studied in detail. These atmospheric 
conditions are simulated and their attenuation characteristics studied. Effect of various 
parameters like wavelength, particle size, distance, and particle size distributions on 
attenuation is also studied. Electromagnetic wave scattering patterns based on the particle 
size and the size parameter are simulated. A simulation package has been developed 
which can calculate the attenuation for a given link under the specified weather 
conditions. We find that the most severe atmospheric conditions are the heavy fog and 
smoke. Heavy rainfall conditions also cause severe losses. 
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CHAPTER 1 
INTRODUCTION 


Over the last decade, there has been a steady increase in the number of 
consumers using high capacity data transmissions, and their data rate demands have 
risen from hundreds of megabytes per second to tens of gigabytes per second. There is 
little room left in the radio frequency spectrum to add more transmitting channel. For 
this reason, many companies turned up towards the light as a way to provide the 
needed room for communication expansion. Also, advancement in optical devices has 
led to an increasing number of applications for optical communication systems. These 
systems consist of both guided and unguided communication links, termed Optical 
Wireless Communication (OWC) Links. 

In contrast to guided systems such as the fiber optic cables used in telephone 
networks, unguided systems launch light into a medium without any attempt to 
control the propagation of the light after it leaves the transmitting optics. Unguided 
systems offer numerous advantages over RF systems including higher data rates and 
less side lobe interference. Also, the narrow beam allows for higher antenna gain and 
smaller antennas. The improvement in system performance over RF designs can be 
used to decrease power, size and weight requirements, an important consideration for 
satellite payloads which are constrained by these parameters. The disadvantages of 
unguided optical systems include the requirement of accurate beam pointing and the 
degradations inflicted by the medium of propagation, such as air or water. Figure 1.1 
is an example of an OWC system. 



Figure 1.1. Example of Optical Wireless Network [ 1 ] 
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The emphasis in this thesis is on the degradations induced by the atmosphere 
on the received optical signal. The major degradations included in this study are: 
1) Attenuation due to fog, 2) Attenuation due to rain, 3) Attenuation due to smog. 
Since outdoor optical hardware and prototypes are expensive and inflexible, 
simulation is an attractive alternative for predicting the performance of optical 
systems. Well established simulation models exist for the guided medium, but for 
unguided systems an end-to-end simulation is very much required. It is the objective 
of this thesis to simulate the various degradations introduced by the atmosphere on the 
performance of outdoor optical wireless links. 

Advantages of OWC systems 

a) Easy to deploy, 

b) no tariffs are required for its utilization (an operating cost saving), 

c) there are no RF radiation hazard (when eye safe power levels are maintained), 

d) has large bandwidth, which facilitates high data rates, 

e) very compact transmitting and receiving modules, possible for short range 
applications, 

f) has low power consumption, 

g) exploits unused part of the electromagnetic spectrum. 

Disadvantages of OWC systems 

a) requires accurate pointing and alignment between transmitter and receiver, 

b) dependent on atmospheric effects - rain, fog, smog and snow outages, 

c) less suited to broadcasting because of the narrow optical beam, 

d) selectivity and availability of sources and detectors in suitable bands. 

1.1 APPLICATIONS OF OPTICAL WIRELESS COMMUNICATIONS 

Optical Wireless communication can bridge the gap by enabling connection with 
offices, business facilities, and other targeted locations without relinquishing the 
performance parameters offered by Radio Communication. Some of the short range, 
long range, and wide area applications are listed below. 
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Short Range Applications 

• Weather monitors; fog, snow, rain using light back-scatter, 

• Traffic counting and monitoring, 

• Optical Interconnects. 

Long Range Applications 

• Deep space probe communications; distances measured in light-years, 

• Building to building very high data rate computer data links; 

• Ship to ship communications with high data rates and security, 

• Telemetry transmitters from remote monitors for weather, geophysical, etc. 

• Electronic distance measurements, such as hand held units out to 1 000 ft., 

• Optical radar; shape, speed, direction and range, 

• Remote telephone links, which are cost effective than microwave links. 

Wide Area Applications 

• Campus wide computer networks, 

• City-wide information broadcasting, 

• Inter-office data links. 

1.2 THESIS OBJECTIVES 

The main objective of the thesis is to study the effect of rain, fog, and smog on the 
performance of outdoor optical wireless links. This involves the following tasks. 

• Review of outdoor optical links. 

• Understanding the various aspects of outdoor optical communication systems 
under different atmospheric conditions. 

• Theoretical study of the drop size distributions under different atmospheric 
conditions. 

• Simulation of various atmospheric losses with different link design 
'parameters. 

• Graphical representation of all the losses studied. 

• Developing a simulation package to calculate the atmospheric losses for the 
given atmospheric condition. 
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1.3 THESIS ORGANIZATION 

•Chapter 2 gives a review of the outdoor optical links. A brief comparison of the 
presently available systems are also discussed. 

Chapter 3 deals with the various design issues and challenges, such as selection of 
wavelength, scattering, absorption, etc., encountered in an outdoor optical wireless 
link. 

Chapter 4 is devoted to the theoretical study of the effect of various atmospheric 
conditions on outdoor optical wireless links, i.e., the effect of rain, fog , and smog. 

Chapter 5 gives the simulation methodologies of outdoor wireless links. The effect of 
the rain, fog and smog on the outdoor optical links is simulated and plotted. An 
analysis on those results is also presented. 

Chapter 6 concludes the thesis with the major results and gives suggestions for 
further work in this area. 
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CHAPTER 2 

REVIEW OF OUTDOOR OPTICAL WIRELESS LINKS 

2.1 INTRODUCTION 

Outdoor optical wireless communication (OWC) systems use light to 
communicate through the air, and require line of sight between the transmitter (Tx) 
and the receiver (Rx) units. Modem systems typically use lasers at the transmitter to 
transmit the signal and photo diodes at the receiver to detect the incoming light and 
send an appropriate signal for processing. Outdoor OWC systems may be used for 
medium to very high data rate applications. Increasingly these systems are used for 
high-bandwidth applications needing data transfer of hundreds of megabits per 
second. Moderately priced OWC consumer products are now available. 

2.2 A BASIC OUTDOOR OPTICAL LINK 

There are three key functional elements of a free space OWC system, viz. the 
transmitter, the atmospheric channel, and the receiver. The transmitter converts the 
electronic signal into light. The light propagates through the atmosphere to the 
receiver, which converts the light back into an electronic signal. The transmitter 
includes a modulator, a laser driver, a LED or laser, and a collimator. The modulator 
converts bits of information into light signal in accordance with the chosen 
modulation method. The driver provides the power for the laser and stabilizes its 
performance; it also neutralizes such effects as temperature and aging of the laser or 
LED. The laser is characterized by its central wavelength, average power, and beam 
divergence angle. Ideally, its frequency spectrum is sufficiently narrow to allow 
optical analysis to relate to the central wavelength alone. The LED is a lower power 
source. It is often preferred to the laser as a cheaper alternative, and is becoming 
increasingly attractive with the recent developments in super-luminescent LEDs with 
relatively narrow bandwidths and up to 70 mW power outputs. The beam divergence 
angle determines the free space power loss, or how large the laser spot will be at a 
given distance from the source. The collimator aligns the laser-LED radiation to a 
collimated beam and directs it to the receiver. 
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The receiver includes a telescope, a filter, a photo-detector, an amplifier, a 
decision device and a clock recovery unit. The telescope collects the incoming 
radiation and focuses it onto the filter. The filter removes background radiation and 
allows only the wavelengths of the signals to pass through it. The photo-detector 
converts the optic radiation into an electronic signal, and the amplifier amplifies the 
electronic signal. The decision unit determines the nature of the bits of information 
based on the time of arrival and the amplitude of the pulse. The clock recovery unit 
operates in parallel to the decision making unit and synchronizes the sampled data 
into the decision making process. The quality of reception is measured by the 
probability of error, expressed in terms of the bit error rate (BER). 


Transmitter Atmospheric Channel Receiver 



Figure 2.1 Block diagram of Optical Wireless Communication Link 

The transmitter and the receiver are separated by the propagation channel, 
which is the atmosphere in the case of OWC. Even on an apparently clear day the 
atmosphere is pervaded by molecules and aerosols, which cause absorption and 
scattering of the light. Changes in the temperature along the propagation path lead to 
scintillations in the received light due to resultant turbulence. 

Line of sight is imperative for outdoor OWC systems, and may be maintained 
with the help of a tracking and pointing system. Tracking is performed in two stages. 
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coarse tracking and fine tracking. Coarse tracking may use GPS or other a priori 
knowledge. Fine pointing requires electro-optic mechanisms such as a quadrature or 
matrix detector. Pointing involves a beam-steering device which may be mechanical, 
such as a galvo-mirror, or non-mechanical, such as acousto-optic crystals, electro- 
optic devices or optical phased arrays. 

2.3 LITERATURE REVIEW 

Free-space optical communication has attracted considerable attention recently 
for a variety of applications. Gagliardi [2] has studied outdoor optical communication 
systems in good detail. Analysis and design of an OWC system is discussed in detail. 
A review of optical fields, sources, channels and signal descriptions is also presented. 

The effects of wind, earthquakes, and thermal expansion on an OWC system 
are described by Kim et al. [3]. They studied the transmission of high speed computer 
network data on the free space channel and also discussed a tracking system to nullify 
the effects of building sway to maintain alignment between the transmitter and the 
receiver. 

Leeb [4] gives the background radiation effects on the transmission of optical 
signal through the atmosphere and studies the degradation of the signal to noise ratio 
with the background illumination. 

Studies based on available weather statistic measurements indicate that the 
weather effects can be mitigated by use of medium- and far-IR laser sources. 
Chapasso et al. [5] discuss the quantum cascade lasers (QCLs), which are adjusted to 
operate within the atmospheric window regions in the mid-wave and long-wave IR 
wavelengths. They demonstrated the use of QC Lasers in effectively mitigating the 
various atmospheric effects. 

'A comprehensive study of the effect of turbulence on communication links is 
described by Brookner [6]. Davis [7] studied the signal to noise ratio for an AM signal 
passing through the turbulent atmosphere and also the probability of detection of 
fluctuating laser signals. 



8 


The effect of fog on the bit-error rate of a free-space laser communication 
system is discussed by Strickland et al. [8]. The paper shows that the BER increases 
as the optical depth (which is defined as the extinction coefficient over the length of 
the channel) increases which in turn leads to increase in the atmospheric attenuation. 

Grotzinger [9] discussed about the performance of free space systems under 
different weather conditions. He also discussed fade margin and gave a review of the 
methods that can be employed to improve the fade margin. 

Kim et al. [10] studied suitable wavelengths for the propagation of laser beam 
in foggy conditions. They showed that in dense fog (i.e. for visibility < 200 m) the 
attenuation of the laser beam is independent of wavelength and all are equally 
attenuated by fog. For light and moderate fog conditions (i.e. visibility < 500 m) 
1550 nm wavelength is less attenuated than the 785 nm wavelength. They formulated 
a new equation for the atmospheric attenuation under different visibility conditions. 

Acampora et al. [11] proposed a hybrid access network that uses small radio 
cells, where the cell base stations are interconnected by optical wireless links. As the 
links are for a short ranges (typically a tens of meters), the system is almost immune 
to all the extreme weather related impairments. 

Achour [12] presented a review of the effects of rainfall on free space optical 
signal propagation. He presented a simulation model for the computation of 
atmospheric attenuation with respect to the rainfall rate. 

Achour [13] presented a review on the effects of fog, haze and low clouds on 
free space optical signal propagation. He presented a simulation model for the 
computation of atmospheric attenuation under the above conditions. 

Satie [14] studied the performance of the outdoor short range (up to 30m) 
optical wireless links under clear weather and fog conditions. He conducted a few 
experiments in fog conditions and established a relationship for the atmospheric loss. 
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2.4 OUTDOOR OPTICAL COMMUNICATION SYSTEMS (PRESENT 

STATUS) 

Optical communications is one of the cornerstones of today’s revolution in 
information technology. With the drive towards portable and multimedia 
communications, we are increasingly faced with the challenge of bringing the 
capacity of our communication infrastructure directly to the user, providing seamless 
access to huge quantities of information, anywhere and anytime. Whether it is the 
transfer of an image from a digital camera to a laptop computer or the communication 
of data within a massively parallel computer, there is an urgent need to develop new 
methods of high speed data communications. Light offers many advantages as a 
medium for communication. Whether traveling through free space or through optical 
fiber, light provides large channel bandwidth and data rates in the range of terabits per 
second. This immense capacity is due to the nature of the photons that constitute an 
optical signal. Unlike electrons, photons react weakly to their environment and to one 
another. As such, optical signals neither generate nor are sensitive to electromagnetic 
interference (EMI), parasitic coupling, and other problems faced by the electrical 
signals. Due to their advantages, optical wireless links are becoming more popular in 
application areas as compared to traditional Radio and fiber-optics. Some of the major 
companies and their outdoor Optical Wireless Systems are listed below. 


• Dominion Lasercom Inc. [15] - 


Product 

Name 

Link 

Distance 

Data-rate 

Wave- 

length 

Type of 

Transmitter 

Output 

Power 

Type of 

Detector 

Skyfiber 

2 km 

1 00 Mbps 

780, 808, 

850nm 

AlGaAs 

Multimode 

30- 

lOOmW 


Skyfiber 

2 km 

1 0 Mbps 

780,808, 

850nm 

AlGaAs 

Multimode 

30- 

lOOmW 



• Laser Wireless Systems [16] 


Product 

Name 

Link 

Distance 

Data-rate 

Wavelength 

Type of 

Transmitter 

Output 

Power 

Type of 

Detector 

ExoLACE 

E155 

3 km 

155 Mbps 

810 nm 

AlGaAs 

Laser Diode 

25 mW 

Si APD 








• Light Pointe Systems [17] 

Product 

Name 

Link 

Distance 

Data-rate 

Wavelength 

Type of 

Transmitter 

Output 

Power 


FlightApex 

1.3 km 

2.5 Gbps 

1 550 nni 

EDFA 


InGaAs 


Flightstrata 


FSA52E, 

FSA52EW 

5.6 -5.2 

km 

1 1 .5 Mbps 

to 54 

Mbps 

FSAI55E, 

1 

OO 

1 .5 Mbps 
to 155 

FSAI55EW 

km 

Mbps 

FSA622 

3.3 km 

622 Mbps 

FSA-G 

2.0 km 

1 .25 Gbps 


850 nm 


VCSEL 


Si APD 


Cable Free Solutions Inc. [18] 




Product Link 
Name Distance 


CF Access 

3 km 

A3 1000 
CF 622 

Series 2 km 

CF 2000 
CF Gigabit 

1 km 

G 1000 


.. , Data-rate Wavelength 

)istance ® 


3 km 155 Mbps 980 nm 

2 km 622 Mbps 780 nm 

1km 1.5 Gbps 780 nm 


Type of Output Type of 

Transmitter Power Detector 


Enhanced 

APD 


Plain Tree Systems [19] 


Product 

Name 

WB500 


WB600 


km 



Data-rate 

Wavelength 

Type of 

Transmitter 

10-100 

Mbps 

800-900 

nm 

LED 

155 Mbps 

800-900 

nm 

LED 
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• MRV Communications Inc. [20] 


Product 

Name 

Link 

Distance 

Data-rate 

Wavelength 

Type of 

Transmitter 

Output 

Power 

Type of 

Detector 

Terescope 

1 55 Series 


TS3I01 


155 Mbps 

850 nm 



Si PIN 

TS 3303 

1.7 km 



TS 4900 

4.2 km 



Terescope 

34 Series 


TS 802/ST 

1.2 km 

34 Mbps 

850 nm 

1 LED 

0.8 mW 

Si PIN 


3 km 

3 VCSEL 

19.5 

mW 

TS 940/ST 

4 km 

4 VCSEL 

24 mW 

Terescope 

10 Series 
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• Fsona Systems Corporation [21] 


Product 

Name 

Link 

Distance 

Data-rate 

Wavelength 

Type of 

Transmitter 

Output 

Power 

type of 

Detector 

8-E 

50 m - 

4.95 km 

1.5-10.5 

Mbps 



100 


155 E 

20 m- 

2.6 km 

1.5-160 

Mbps 



mW 


155-E 

100 m- 

3.8 km 

31-180 

Mbps 



320 

mW 


1250 E 

100 m - 

3.2 km 

100- 

1448 

Mbps 

1550 nm 

Laser Diode 

280 

mW 

- 

2 M 

300 m- 

7.7 km 

31 -62 

Mbps 



640 

mW 


1550 M 

400 m- 

5.3 km 

100- 

1448 

Mbps 



560 

mW 



• Laserbit Wireless Cables [22] 


Product 

Name 

Link 

Distance 

Data-rate 

Wavelength 

Type of 

Transmitter 

Output 

Power 

Type of 

Detector 

LB 1500 

1.5 km 




2x70 

mW 


LB 2500 

2.5 km 

155 Mbps 

785 nm 

Laser Diode 

4x70 

mW 

APD 

LB 5000 

5 km 




8x70 

mW 



The commercially available links are broadly divided into two categories. 
Short range links (having link range less than 2 km) and long range links (having link 
range greater than 2 km). 

Short Range Links 

Light Pointe Communications Inc. has introduced two short range products 
FlightApex and FlightStrata for up to 1 .3 km and 2 km link range and 2.5 and 1 .25 
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Gbps data rates, respectively. Cablefree Solutions Inc. has also launched a product 
range CF Gigabit, having a link range of 1 km operating at a data rate of 1 .5 Gbps. 
Canonbeam from Canon and CF 622 product range from Cable free Solutions Inc. is 
available for link ranges up to 2 km at 622 Mbps. For the speed of 155 Mbps, many 
companies viz. Cablefree Solutions, Fsona Systems Corporation, Laserbit Wireless 
Cables and Dominion Inc. have a wireless link up to a distance of 2 km. While most 
products are only for point to point communications, companies such as AirFiber and 
Terabeam have brought out products that easily allow a mesh of links to be set up. 

Long Range Links 

There are a number of long range products available for the different data rates 
for link ranges up to 7 km. Fsona System Corporation’s M and S Series range of 
products provide links from 3.3 km to 7.7 km range for data rates starting from 
31 Mbps to 1.448 Gbps. For the link range of 3.5 km, LightPointe 
Communications Inc. provide data rate of 622 Mbps. For the speed of 155 Mbps and 
lower, many companies viz. Cablefree Solutions, Fsona Systems Corporation, 
Laserbit Wireless Cables and Dominion Inc. have wireless links up to a distance of 
5km. The power level employed for the long distance links generally fall into class 3 
eye-safe. LSA Photonics has introduced a product for relatively lower data rate (45 
Mbps) but longer ranges (up to 15 km). 

In this chapter a brief review of all the advancements that has taken place in 
this field is given. A comprehensive review of all the existing outdoor links has also 
been presented and compared. In the next chapter, the design challenges that are faced 
by the OWC system are discussed. 
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CHAPTER 3 

DESIGN CHALLENGES IN 
OUTDOOR OPTICAL WIRELESS LINKS 

With the optical wireless communication becoming very popular, there is a 
need to build robust optical transmitters and receivers for establishing point to point 
communication through the atmosphere. Despite the advantages, there are some 
important design issues associated with the OWC systems. Some of the major design 
challenges faced by an OWC system are selection of wavelength, selection of source 
and detector, problem of alignment between transmitter and receiver, atmospheric 
effects such as absorption, scattering and turbulence. These issues are briefly 
discussed in the following section. 

3.1 SELECTION OF WAVELENGTH 

Performance of OWC systems depends a lot on the wavelength used for 
communication, thus making the selection of wavelength a crucial issue. The two 
factors which needs consideration while specifying the operating wavelength are the 
eye safety issue and the atmospheric attenuation. 

The system must not pose any danger to the people who encounter the 
beam. Therefore a laser safety standard has been established in which optical sources 
have been classified in accordance with their total emitted power. The principal 
classifications are summarized in Table 2.1 for a point source emitter such as the 
semiconductor laser [23]. 

Outdoor point-to-point systems generally use high power lasers that operate in 
the class 3B band to achieve a good power budget. The 780 to 850nm wavelength 
band is presently the best choice for most applications of optical wireless links, due to 
the availability of low cost laser diodes (LDs). It also coincides with the peak 
responsivity of inexpensive, low capacitance silicon photodiodes. The primary 
drawback of radiation in this band relates to the eye safety; as the radiation at this 
band can pass through the human cornea and be focused by the lens on to the retina, 
causing thermal damage. The cornea is opaque to radiations beyond 1400nm, 
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considerably reducing potential ocular hazards, so that it has been suggested that the 
1550 nm band may be better suited for infrared links. It may be noted that the 
allowable safe power at 1 550nm is about 50 times higher than that at 800nm [1]. 


Wavelength 

650nm 

SSOnm 

ISlOnm 

1550nm 

Classl 

Up to 0.2mW 

Up to 0.5mW 

Up to 8.8mW 

Up to lOmW 

Class2 

0.2-1 mW 

N/A 

N/A 

N/A 

ClassSA 

l-5mW 

0.5-2.5mW 

8.8-45mW 

10-50mW 

ClassSB 

5-500mW 

2.5-500mW 

45-500mW 

50-500mW 

Table 3.1 Laser safety c 

assifications for a point source emitter 


Second issue is the atmospheric attenuation. Atmospheric optical 
attenuation processes are wavelength dependent, such as Rayleigh scattering from air 
molecules which varies as the fourth power of wavelength. The other attenuation 
process is Mie-scattering which causes more attenuation than Rayleigh scattering. The 
wavelength dependence of Mie scattering is sensitive to the nature of fog droplets; 
stronger when the size of the droplet is nearly equal to the wavelength of light. But in 
light fog and haze the attenuation due to Mie scattering is less. The wavelength 
dependence of the scattering coefficient P is given by [1 1]. 

3.91 r ^ 


P = - 


550nm 


(3.1) 


where V is visibility in km, and q = 0.585 , for V<6 km. Hence for OWC links used 

in high fog situations, lower wavelength source will be the right choice. 


Another factor affecting the wavelength selection is the beam spreading due 
to diffraction which is linearly proportional to the wavelength [23]. The resulting 
attenuation can be given by (4 ;t A, / where A, is the transmitting area. 


For terrestrial OWC links the other challenging issue is the background 
radiation due to sunlight [23]. The solar background radiation is four times less 
effective for 1 550nm than for 800nm .The performance of receiver is also wavelength 
dependent. Generally the photodiodes near 800 nm and 1550 nm achieve comparable 
quantum efficiency. The output response of the detector depends on the energy of the 
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photon. A photon at SOOnm has roughly double the energy of that at 1550 nm. Hence 
for the same noise conditions, the optical energy at 1550 nm will be approximately 
3dB less. This necessitates higher receiver sensitivity in the case of 1 550 nm systems. 
Also the capacitance and gain of the APD are also wavelength dependent quantities. 
Hence the selection of wavelength for an OWC system should consider all above 
factors. 


3.2 SELECTION OF SOURCE 

Typical lasers used in OWC are edge-emitting semiconductor devices, with 
wavelengths 785, 850, 980 or 1550 nm. Various output powers and wavelengths are 
available. Lasers have the advantages of high output powers, high modulation rates, 
and if cooled, long-term reliability. Nowadays, VCSELs (vertical-cavity surface- 
emitting lasers) are also becoming very popular as transmitting devices [18]. 



p AIGaAs 


nAIGaAs 


n GaAs Substrate 


Top Metaffic 
Contact 


Proton Bombarded 
Semi-lnsulafing Barrier 
p+ GaAs 


Bottom Contact 


Emitting Region 


Light Out 



Au/Zn Circular 
mirror and electrode 


Active 

Region 

Si02 


Edge Emitting Laser VCSEL 

Figure 3.1 Types of sources for OWC systems 

VCSELs are designed using a vertical cavity, and as a surface emitter, emits light 
out of its top surface. This sets VCSELs apart from all other lasers, classed as edge- 
emitters, which have a cavity that lies in the plane of the wafer (usually thought of as 
horizontal) and emits light out of the side. 


VCSELs operating at 850 nanometers have become popular in Local Area 
Networks due to the following factors. 

• At 850nm, VCSELs are cheap, and offer high bandwidths. 
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' • VCSELs emit a uniform, narrow, circular beam, which simplifies optical 

system design considerably. In contrast, edge-emitters have an elliptical beam, 
which is strongly divergent. 

• A VCSEL laser cavity is very short, compared to that of an edge-emitter. So, 
when the laser is turned off, it takes less time for the remaining light to exit the 
cavity. As a result it is possible to modulate the laser (turn the beam on and 
off) directly, which avoids the need for an external modulator. Devices w'ith 
direct modulation speeds of 2.5 Gbps are common place. 

• VCSELs are very efficient at turning injected carriers into photons. Typical 
efficiencies are in the region of 6 to 25 per cent. As a result they consume less 
electrical power than an equivalent DFB (distributed feedback) laser. For 
example, a DFB laser might draw 60 mA, while a VCSEL would only require 
15 mA to produce the same optical output. 

The problems associated with the use of VCSELs are: 

• VCSELs have an abrupt failure mechanism. Long-term-reliability at 
moderate-to-high temperatures is a problem, and hence VCSELs have been 
dominant in LAN environments and not in telecoms, where high reliability is 
paramount. 

• The materials used to make 850nm VCSELs are not yet applicable to longer 
wavelengths. Currently, commercial VCSELs that operate at the wavelengths 
of 1 3 1 0 or 1 550 nm have yet to become generally available. 

LEDs are also used in some OWC links but only for very low cost short range 
systems. 

3.3 SELECTION OF DETECTOR 

There are a number of factors to be considered while examining the 
effectiveness of the receiver in an OWC system; these include the type of detector 
used, the sensitivity rating and size of the detector, the size and design of the receiver 
optics [18]. Optical detectors used in OWC equipment come in two basic types: PIN 
and APD. The PIN detector is a lower cost detector that has no internal gain, while the 
APD is a more expensive but also more sensitive detector with internal gain. 
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Parameter 

PEV Photodiodes 

APDs 

Construction Materials 

Si, Ge, InGaAs 

Si, Ge, InGaAs 

Wavelength 

0.6 to 1 .8 pm 

0.6 to 1.8 pm 

Responsivity 

0.5 to 0.85 Amps/Watt 

0.5 to 100 Amps/Watt 

Gain 

1 

10-150 

Support Circuitry Required 

None 

High Voltage, 

Temperature Stabilization 

Cost 

Low 

High 


Table 3.2 Comparison of PIN photodiode and APD 


One of the key parameters to consider when selecting a detector is spectral 
noise. A detector normally operates in one of two noise-limited detection regimes; 
either detector noise limited at low power levels, or photon shot noise limited at 
higher power level. In an APD, sensitivity at low light levels will be limited by the 
shot noise and the APDs leakage current. Shot noise derives from the random 
statistical fluctuations of the dark current, and signal current. Dark current or shot 
noise in a PIN differs for an APD, as in an APD the bulk leakage current is multiplied 
by the gain, M, of the APD. In addition, the avalanche process statistics generate 
current fluctuations, and APD performance is degraded by an excess noise factor (F) 
compared to a PIN. At higher signal light levels, the detector transitions to the photon 
shot noise limited regime. In the absence of other noise sources, an APD provides a 
signal-to-noise ratio (SNR) which is worse than a PIN detector with the same 
quantum efficiency. An APD, however, can produce a better overall system signal-to- 
noise ratio than a PIN detector in cases where the APD internal gain boosts the signal 
level without dramatically affecting the overall system noise. In summary, an optical 
receiver using an APD offers superior sensitivity than PIN at medium-to-high 
bandwidths. A choice of Silicon, Germanium or InGaAs offers operation at a variety 
of wavelengths. 

Another important consideration is the size of the detector or the active area. 
Active area corresponds to the size of the actual light sensitive region, independent of 
the package size. PIN Diodes with large active areas will capture more light but will 
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always be slower than smaller devices and will also produce more noise. However, if 
a small device contains an attached lens it will often collect as much light as a much 
larger device without a lens. However, the devices with attached lenses will collect 
light over narrower incident angles (acceptance angle). Flat surface devices are 
usually used if light must be detected over a wide area. For most applications either 
style will work. For high speed applications a device with a small active area is 
always preferred. However, there is a tradeoff between device speed and the active 
area. For most long-range applications, where a large light-collecting lens is needed, a 
large area device should be used to keep the acceptance angle from being too small. 
Smaller acceptance angles can make it nearly impossible to point the receiver in the 
right direction to collect the light from the distant transmitter. 

3.4 ALIGNMENT PROBLEM 

Continuous alignment between transmitter and receiver during communication 
is essential for the successful data transmission. This is particularly problematic when 
narrow beam divergence angles and fields of view (FOVs) are used. Thermal 
expansion of building frame parts and weak earthquakes are the possible causes of 
misalignment [1]. Another common cause of misalignment is wind, especially when 
communication equipment is placed on tall buildings. Under the influence of dynamic 
wind loads, many high-rise buildings sway in along-wind and across-wind directions, 
or twist due to torsion. This interrupts the line-of-sight. The horizontal movement of 
the building due to these effects can vary from 1/200 to 1/800 of the building height. 
Building sway is a random process that affects system performance since it is a source 
of pointing error at the beam steering stage. Knowledge of the stochastic behavior of 
the pointing error enables corrective transceiver design or solutions on the system 
level such as intelligent protocols. 

3.5 SCATTERING 

It is generally agreed that the primary drawback of OWC system is the 
downtime that results from adverse weather conditions. Although rain and snow can 
hamper light propagation, the most detrimental environmental conditions are fog and 
haze. This is due to their high concentration of scattering particles with radii of the 
same order of magnitude as the laser wavelength [1]. Scattering is the deflection of 
incident light from its initial direction, which causes spatial, angular, and temporal 
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spread. The precise scattering mechanism is dependent on the ratio of particle size to 
radiation wavelength. When this ratio is near to unity, the direction of propagation 
subsequent to scatter is determined by the Mie-scattering solution to the Maxwell 
wave equation. Given the size distribution of the scattering particles, the wavelengths 
of radiation, and the refractive indices of both the drop and the air, the probability 
function of the Mie-scattering angle can be derived. Hence, the scattering and 
resultant communication degradation can be predicted. Raindrops and snowflakes are 
very much larger than optical wavelengths, and are not subject to Mie-scattering and 
do not seriously impede light propagation. Scattering results in attenuation and 
distortion of the transmitted signal. Attenuation of the order of tens or even hundreds 
of decibels per kilometer can be encountered when the propagation channel is 
unfavorable. 

3.6 ABSORPTION 

Molecules in the atmosphere absorb laser light energy, causing attenuation of 
the optical power [1]. This phenomenon is extremely sensitive to wavelength so that a 
given environment may be transparent to some wavelengths, while almost entirely 
blocking others. The infrared absorption by the atmosphere is largely due to water 
vapor and CO 2 , while below around 200 nm absorption by O 2 and O3 are so high that 
there is virtually no transmission. Regions of fair transmission are called free 
absorption windows, and are found in the visible and near infrared range. At middle 
infrared wavelengths there are more absorption windows in the region of 3-5 pm and 
8-12 pm, but above 22 pm the water vapor absorption again rises to prohibitive levels 
for transmission. Weather can often change rapidly, and as the molecular content of 
the atmosphere alters, the transmission characteristics of the propagation channel 
change too. These changes are due to absorption as well as scattering and may 
seriously influence the availability of a link. Paradoxically, the absorption of the 
transmitted signal can be seen as an advantage when strict security is a priority, as 
beyond-a limited range no power can be detected. 

3.8 TURBULENCE 

Turbulence effects are due to varying temperatures and humidity levels within 
the propagation media. It causes scintillation of the light and “beam wander” [6]. In 
other words, the amplitude of the received signal fluctuates over time, and the 
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position of the incident light shifts in space. This may be due to air conditioning vents 
in the vicinity of the transceiver, irradiative heat from the roof or currents of 
pollutants in the atmosphere. Temperature and humidity gradients cause changes in 
the atmospheric refractive index, which is the source of optical distortions. Winds and 
cloud coverage also influence the level of turbulence, and even the time of day can 
alter temperature gradients, and hence turbulence effects. A number of statistical 
models have been developed to describe the dependence of the refractive index on the 
ambient temperature, pressure, and humidity as well as on the radiation wavelength. 
The basic statistical parameter commonly used for analyzing power transmission is 
the spatial cross-correlation of the randomly varying refractive index. The power 
spectral density of the index fluctuations defined for spatial frequency is the three- 
dimensional Fourier transform of the spatial cross-correlation function. To model 
turbulence, different energy exchange mechanisms have been distinguished, 
corresponding to different ranges of spatial frequency. Commonly, the strength of the 
turbulence is described using the auxiliary term (the refractive index structure 
constant) where inertial forces dominate. This inertial range of spatial frequencies is 
called the Kolmogorov spectrum and has been extensively investigated both 
theoretically and empirically. Using Ritov theory, a lognormal power intensity 
distribution function has been derived for the spatial distribution of light propagation 
through a turbulent atmosphere. 

In this chapter, the major design challenges in outdoor optical wireless link 
have been briefly discussed. In the next chapter, the theory of radiative transfer and 
scattering is discussed in detail. Also, the various distributions available for rain, fog 
and smog are discussed. 
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CHAPTER 4 

THEORETICAL STUDY OF OUTDOOR OPTICAL WIRELESS 
LINKS IN RAIN, FOG AND SMOG 

4.1 THE ATMOSPHERIC CHANNEL 

The atmosphere is composed of collections of gases, atoms, water vapor, 
pollutants and other chemical particulates that are trapped by the Earth’s gravity field and 
it extends to approximately 400 miles in altitude. The heaviest concentration of these 
particles is near Earth in the troposphere level with particle density decreasing with 
altitude up through the ionosphere. Actual particle distribution depends on the 
atmospheric conditions. These atmospheric particles interact with all radiation fields that 
propagate through the radiation belts, with the primary effects being power losses and 
wavefront distortion. 

Power losses and distortion are caused by the absorption and scattering of the 
radiation fields by the particulates. These effects become most severe as the radiation 
wavelengths begin to approach the cross-sectional dimension of the particulates. Because 
particle sizes in the atmosphere range from centimeters down to micrometers, the 
atmosphere is especially deleterious to optical transmissions. 

Atmospheric conditions can be roughly classified into three basic types: clear air, 
clouds and rain. The clear air channel is the most benign, characterized by long range 
visibility, clear weather and relatively low attenuation. Clear air, however, can still 
contain eddies and temperature gradients (clear-air turbulence) that can produce changes 
in the index of refraction of impinging fields. These index changes act as optical lenses 
that refocus and reorient beam propagation. Cloudy atmospheric conditions can range 
from mist and fog to heavy cloud cover extending from near Earth to upper altitudes and 
is characterized with increasing water vapor accumulation and higher levels of 
attenuation. Rain represents the presence of water droplets of significant sizes, and can 
produce the severe effects, depending on rainfall rate and rain cloud extent. 

As field radiation impinges on an atmospheric particulate, a portion of its energy 
is absorbed, and the angle of the remainder is redirected. The particle absorption produces 
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a field power loss and the angle redistribution produces scattered energy. Surrounding 
particles further absorb and scatter the redistributed field, producing aggregate field 
attenuation and scattering mechanism on the overall field as it propagates. It is precisely 
this scattering within a sunlit cloud, for example, that makes it appear uniformly 
illuminated when observed from below. 



0.6 10 1.5 2.0 2.S 3.0 30 4.0 4.5 5.0 5.5 


Figure 4.1 Transmittivity of the atmosphere [2] 

Figure 4.1. shows a plot of the field transmittance observed at various optical 
wavelengths propagating through the clear air atmosphere. The significant attenuation 
observed at selective wavelengths indicates field absorption from specific particles and 
shows the importance of properly selecting wavelengths in the high transmittive bands for 
space optics systems. 

The important parameters that represent the atmospheric conditions are the 
particle size (cross sectional dimension relative to the wavelength), and the particle 
density (volumetric concentration of particles). Low densities of large particles tend to 
obey single scatter theory and produce primarily attenuation with mostly forward 
scattering at angles close to the normal. High densities of small particles tend to obey 
multiple scatter theory, which produce elemental beam scattering and random phase 
shifting across the wave-front as it propagates. 

Figure 4.2 shows typical values for average particle size and the corresponding 
particle densities for various cloud and rain conditions. It is evident from the figure that 
the conditions vary from high density small particles during mist and fog to low-density 
large-particle concentrations during heavy rain. 
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ipm) (cm) 

Droplet radius 

Figure 4.2 Cloud and precipitation drop size distribution [2] 

The different types of atmospheric constituent’s sizes and concentrations can be tabulated 



Table 4.1 Atmospheric constituents, their sizes and concentrations [12] 
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4.2 ATMOSPHERIC CONDITIONS 

Under different atmospheric conditions, an optical wireless communication system 
undergoes different attenuation in accordance with the theory of radiative transfer (to be 
discussed in the following section). The attenuation can vary from 0.4 dB/km to 500 
dB/km, depending upon the weather conditions. The laser beam propagating through the 
atmosphere is attenuated by scattering and absorption of radiation by the following 
atmospheric conditions: 

1. Rain 

2. Fog 

3. Smog 

4. Dust 

The impact of the various weather conditions on an OWC link is shown in Fig 4.3. It can 
be seen that in heavy fog, the attenuation over a link can be a few hundreds of dB/km, 
whereas that for rain it is a few tens of dB/km. 
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Figure 4.3 Impact of weather on laser power 


4.3 SCATTERING 

When a beam of light impinges on a particle, electric charges in the particles are 
excited into an oscillator}' motion. The excited electric charges reradiate energy in all 
directions (scattering) and may convert a part of the incident radiation into thermal energy 
(absorption) [I6]and[17]. 
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The energy scattered by a particle (Fscat )is proportional to the incident intensity F, and is 
given by, 

F scat ~ CscatFi (4.1) 

where Cscai has units of m\ is the single-particle scattering cross section. For absorption, 
the energy absorbed (Fabs ) is described analogously as, 

P abs ~ CabsFi (4.2) 

where Cabs in units of , is the single-particle absorption cross section. 

Conservation of the energy requires that the light removed from the incident beam 
by the particle is accounted for by scattering in all directions and absorption in the 
particle. The combined effect of scattering and absorption is referred to as extinction, and 
an extinction cross section Cext can be defined by 

Cext ~ Cscat Cabs (4.3) 

The dimensionless scattering efficiency of a particle Q^,, is CscaM, where A is the cross- 
sectional area of the particle. Defining i^ofeand Qext in the same way, 

Qext ~ Qscat Qabs (4.4) 


Incident Beam wth wavclcngtli X Incident Beam with wavelength 



(c) Small Particles with diameter > X 


Figure 4.4 Scattering patterns of electromagnetic waves by spherical drops 
(a) Rayleigh Scattering (b) Mie Scattering (c) Geometric Scattering [24] 



Light scattering mechanisms of the particle can be divided into three categories 
depending upon the size parameter x (which is equal to 2itr/X, wiiere r is the radius of the 
particle) and the wavelength : 

1. Rayleigh Scattering 

2. Mie Scattering 

3. Geometric Scattering 


43.1 Rayld^ Scattering 

When the particle is much smaller than the wavelength of the incidait ligiit, 
Rayleigh scattering takes place [25], With respect to light in the visible range of the 
spectmm, particles of diameter <0.1 pm lie in the Rayleigh scattering regime. In this 
regime the pattern of scattered light intensity is symmetrical in the forward and backward 
directions and more or less indep«ident of particle sh^e. Rayleigh scattering occurs 
primarily off the gaseous molecules in the atmosphere. The attenuation coefficient varies 
as ( where X is the wavelength). 

Incident li^t on a particle having very small size compared to the wavelaigth, 
behaves in the same way as if the particle is placed in a homogeneous electric field. The 
particle’s own field, caused by the electric polarization of the particle, modifies this field 
both inside and near the particle. This induces a dipole moment. The oscillating dipole 
radiates in all directions. This type of scattering is called Rayleigh scattering. 


Two conditions required for Rayleigh scattering are: 

1 . Size of the particle should be small compared to the wavelength 

2. The applied field should penetrate so fast into the particle that the static 
polarization is established in a time short con:q)ared to the period. 


For an incident wave of unit density, the intensity of the scattered wave at a distance / 
from the particle is given by [25] 


_ 16;7V* f e, 

^ + 2 ^ 2 ) 


sm y/ 


(4.5) 


or 
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where m is the refractive index, X is the wavelength in the medium and ^ is the angle 
measured from the scattering direction to the dipole. 

The extinction and scattering efficiencies for particles in the Rayleigh scattering 
regime are given by [25] 
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and 


QscaXfn,x)=^X* 


- 1 


m' +2 


(4.8) 


4.3.2 Mie Scattering 

As the particle size approaches the laser wavelength, the scattering of radiation off 
the larger particle becomes more dominant in the forward direction as opposed to the 
backward direction [25]. This type of scattering, where the size parameter x varies 
between 0.1 and 50 is called the Mle scattering. The lasercom wavelengths are Mie 
scattered by haze and smaller fog particles. The detailed explanation of Mie Scattering is 
given in section 4.6. 

4.33 Geometric Scattering 

When the particle size parameter x is greater than 50, the scattering is called 
geometric or non-selective scattering [25]. This type of scattering is strongly dependent 
on the particle shape and orientation relative to the incoming beam. The scattering is 
called non-selective because there is no dependence of the attenuation coefficient on laser 
wavelength. 

The fundamental fact is that the incident beam of light, which forms a plane 
wavefront of infinite extent, may be thought to consist of separate rays of light which 
pursue their own path. A small part of a wide wavefront can be looked upon as defining a 
ray which has for some length its own existence independent of the entire wavefront. A 
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length / requires a width of the order of V(/A,), and the requirement for any separate 
existence at all is a width > X. For a particle with a size of 20 or more times the 
wavelength, it is possible to distinguish fairly sharply between the rays incident on the 
particle and the rays passing along the particle. Among the former it is possible to 
distinguish rays hitting various parts of the particle’s surface. 

The rays hitting the particle and passing along it give rise to two distinct phenomena: 

1 . Reflection and Refraction 

Rays hitting the particle’s surface are partially reflected and partially refracted. 
Refracted light may emerge after another refraction, or after several internal 
reflections. The light so emerging and the light directly reflected from the outer 
surface both contribute to the total scattering by the particle. The energy that does 
not emerge is lost by absorption inside the particle. 

2. Diffraction 

The rays passing along the particle forms a plane wavefront from which a part, in 
the form and size of the geometrical shadow of the particle, is missing. This 
incomplete wavefront gives rise to a certain angular distribution of intensity, 
called the Fraunhoffer diffraction pattern. 

4 
3 

O 2 
1 
0 

0 10 20 30 40 

r = radius {j-im) 

Figure 4.5 Scattering efficiency of water particles scattering 
785 nm light as a function of particle size 

As the particle size increases, the value of Qext approaches to 2. The reason for this 
is as follows. The efficiency represents the ratio of energy removed from the wave to the 
total energy intercepted by it. All energy falling on the particle is scattered or absorbed, 
i.e., in any event removed from the proceeding wave. This gives an effective cross- 





30 


section equal to G, the geometric area. Besides that we have the diffraction forming an 
angular pattern that is identical with the diffraction through a hole of area G by Babinet’s 
principle. This gives diffraction with a cross-section G. The total energy removed from 
the proceeding wave corresponds to a cross-section 2G, i.e., to an efficiency factor 

Qext=2. 

4.4 RADIATIVE TRANSFER 

The basis of classical transfer theory is the Lambart Bouger law, which states that the 
change in intensity of radiation due to an interaction with matter depends directly on both 
the incident intensity and the amount of matter [26]. A beam or pencil of radiation is 
attenuated by its interaction with the medium through which it passes. If the initial 
intensity is 1, then after traveling distance ds along the direction of propagation, the 
change in the intensity will be 

dl = {s)plds = -P,, {s)Ids (4.9) 

The parameter p is the density ig/cm^) of the attenuating medium. Equation (4.9) defines 
the mass extinction coefficient kgxt and the volume extinction coefficient Pext. Pext has the 
dimensions of cm"', and kext have the dimensions cm^/g. 

The theory also assumes that emission by an element of mass or volume is linear 
in both intensity and the amount of matter present. The net increase in intensity in a 
distance ds along the direction of propagation is s(s)ds. The equation of transfer for a 
given path through a medium can be thus expressed as 

^ = (5)7(5) + ^(5) (4.10) 

ds 

Formal integration of this one-dimensional transfer equation of the path from s = 0 to s = 
S gives 

.S’ 

7(5) = 7(0)e'^ + jff(5')e-<^-^'’£/5' (4.1 1) 

0 

in which th^ optical thickness along the path is r 

.V 

= \l3exXs)ds 

0 
s 

0 


(4.12) 
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This formal solution to the equation shows that the resultant intensity is the sum of two 
distinct terms; the incident intensity attenuated by extinction along the path (0,S) and a 
component due to emission at points along the path, attenuated along the path (s ,S). If we 
neglect the emission term, we have the usual form of Lambart-Bouger law. 

/(/?) = /(0)exp(-/?„,/?) (4.13) 

where pext is the mean value over the path and R is the distance over which the intensity 
reduction is to be calculated. 

Beer— Lambert Law 

The alternative form of the radiative transfer theory is the Beer-Lambart law also 
known as Beer’s law or the Beer-Lambart-Bouger law. This law is an empirical 
relationship in relating the extinction of light to the properties of the material the light is 
traveling through. 

It states that the extinction coefficient p(X), of a species at a particular wavelength 
of electromagnetic radiation X, is proportional to the concentration c of the scattering or 
absorbing species, of the electromagnetic radiation through the sample containing the 
absorbing species. This can be written in the form: 

p{X) = a{X)c (4.14) 

The proportionality constant a(X) is called the extinction cross section of the 
species at the wavelength X. 

The value of the extinction coefficient p varies between different absorbing and 
scattering materials and also with wavelength for a particular material. It is usually 
determined by experiment. 

The law tends to break down at veiy high concentrations, especially if the material 
is highly scattering. If the light is especially Intense, nonlinear optical processes can also 
cause variances. 

Extinction coefficient is a measure of the amount of radiated power absorbed and 
scattered. Mathematically, 


P{X):=\Og{PJP) 
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where Po is the incident power and P is received power. 
Hence, Beer law alternatively can be stated as 

PiR) = P(0)exp(-o-cl) 

P(i?) = P(0)exp(-AJ) 


4.5 VISIBILITY 


This is an important parameter needed to calculate the optical attenuation due to 
Fog, Haze and clouds. Visibility is generally expressed as a distance. It is defined as the 
maximum distance at which the object contrast against its background is 2% using 550nm 
wavelength light that the human eye is most sensitive to. If the contrast at a distance z is 


C(z), then the relative contrast, defined as the ratio of C(z)/C(0) can be expressed 


C(z) _ 


4.„(0) ^ 

C(0) 

^ ^min (^) / 

^ ^max (0) - -^min (0) > 


as [27]. 
(4.16) 


where Lmax(z) refers to the object luminance and is background luminance at 

C(z) 

distance z. The value of z for = 0.02 is called the visual range (V). If 

f'/'ot * V /• o we assume 


the background irradiance to be zero at z=0 and z=z, then Eqn.(3.28) reduces to 

C(z) L^^jz) 

C(0) L_(0) 


From Ref [27] we see that Lmax(z) can be approximated as 

L^.A^) = L^^iO)cxp{-/3,,,s) 
where a is the attenuation coefficient. 

Hence at visual range (i.e. at s=V) we get 

= exp(-P,,,V) = 0.02 

C(0) 

Hence 

/?^.^,F = -ln(0.02) = 3.912 
Substituting, we get Fin kilometers as 

|ln(0.02)| 3.912 

ft„(^=550mi) /S„{A = 550nm) W'I7) 

If the visibility for an outdoor optical link is given, then the extinction coefficient can be 
calculated using the following equation 


f 


a 


p= 


3.91 

V 


v» 


0.55 /jm 


isiiere q = 1.6 for good visibility i.e. V > 50 km, 
q = 1.3 for 6 km < V < 50 km, 
q = 0.585 f^^^forV<6km 


(4.18) 



Figure 4.6 Attenuation versus visibility [13] 

If we look at Fig 4.4, we will find that dK attenuation value based (m equation 
4.18 is higher than the value calculated by using the drop distribution model It is due to 
the independence of, exponent value and the one-to-one rdation between visibility and 
attenuation coefficirait, firom droplet rizes and disttibutions. As die drop rize distribution 
varies firom place to place, the value of attamadon can not be same. 

To calculate the attenuation coefficient accurately, one can make use of the Mie- 
Debye theory (discussed in the next section). 

4j6 MIE-DEBYE TBKORY 

Mie Debye theory is the rigorous treatmoit of electromagnetic wave scattering by 
s^dierictd particles in a surrounding medium [28], 
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Assumptions made while using Mie Debye theory: 
1 . Particle is a sphere. 


2. Particle is homogeneous (therefore it is characterized by a single refractive index 

m = n-ik at a given wavelength). 


The size parameter in Mie theory is defined by x= InrlX where r is the radius of drop and 

X is the wavelength of incident light in the medium. Relative refractive index of the 

sphere is m = risphere / n,„edium with risphere the complex refractive index of the sphere and 

nmedivm that of the medium. Mie coefficients are expressed as 


^ {n'‘x)if/„ (x) - m {Tnx)y/n (x) 

il/„ {mx)i^ n (x) {mx)^ n (x) 

(4.19) 

^ m - ^„(7Mx)'P’„(x) 

m If/ „ (mx)C„ (x) - if/„ {mx)C'n (x) 

(4.20) 

where if^„ix) = xj„(x) and 4;(x) = x/z^^^«(x) = x( 7 ;(x)-/«„(x)) are the 

functions with the following properties, 

Ricatti-Bessel 

W = (2« + Wnix)/x - lf/„_^ (x) 

(4.21) 

C„,,(x) = (2n + \)C„(x)fx-C„.,ix) 

(4.22) 

If/' „ (x) = if/„,^ (x) -rnf/„{x) lx 

(4.23) 

^'"ix) = Cn-^i.x) - n(;„{x)lx 

(4.24) 

^i/'_l(x) = C0s(x) 

(4.25) 

^o(x) = sin(x) 

(4.26) 

C-i (x) = cos(x) + i sin(x) 

(4.27) 

Co (^) = cos(x) - i sin(x) 

(4.28) 

The index runs from « = 1 to infinity, but the infinite series in Mie formulas can be 

truncated to a maximum rimax, as proposed by Bohren and Huffman, 


nmax = 2 + X + 4x'^^ 


The scattering amplitude functions are 


‘^i(®) = S . ,Ja«^„(cos0) + 6„r„(cos0)] 

(4.29) 

‘S' 2 (®) = Z , ,j6„a-„(cos0) + a„r,.(cos©)] 

^77(n + l) 

(4.30) 


where 7t„ and r„ are Mie angular coefficients. 
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n (cos0) = -: P'„(cos©) 

sin0 


(4.31) 


r (cos0) = — P'„(cos0) 
oQ 

where P\, are the associated Legendre Polynomial. 

From equation 4.26, 

P'„(cos0) = = sin0;r„(cos0) 

a© 


and by differentiation 

dP„{cos<d) 


= cos©;r„(cos©)-sin^ ©;rn(cos0) = r„(cos0) 


The recurrence relations of Legendre polynomials are 


P„(cos0) = COS0 


r2w— ii _ . 

P„_i(cos©)- 

_ n — \ _ 


P„_2(cos0) 


(4.32) 


P „ (cos 0) = (2« - 1)P„_, (cos 0) + P „_2 (cos 0) 

For the far field zone, the solution of the vector wave equation can be obtained as 

E'/ _ exp(-ikr + ikz) 5'2(0)O E'l 
E'r ikr 0 S',(0) p', 


(4.33) 


(4.34) 


(4.35) 


(4.36) 


where £'/ and E'r are the parallel and perpendicular components of the incident electric 
field, and Efi and Efr are the parallel and perpendicular components of the scattered 
electric field. 


The scattered far field for a unit amplitude incident field (where the time variation 
exp(icot) has been omitted)is given by 

P'/ = cos{(t>)S., (cos©) (4.37) 

ikr 

E% = ^^^ -^ ~^^^ -sin(^)P,(cos0) (4.38) 

-ikr 

Intensity for the natural light of unit intensity 

I = («») 

Stt /- 

where /, = |iS',(0)'|and /, = |■S'2(0)^| are the intensity functions. 

Mie efficiency factors are derived from the Mie scattering amplitude 
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a., =4S(2« + l)ReK+A„] (4.40) 

^ /t=l 

Q.a.=\Yi'^ri + \)[\aJ +\b„f] (4.41) 

^ «=1 

a/,,=a..-ac,„ (4.42) 

where Qexi, Qscai and Qabs are the extinction efficiency, scattering efficiency and 
absorption efficiency respectively. As x — > oo, the extinction efficiency Qexi approaches to 
a constant value 2. 

The coefficients of extinction, scattering and absorption can be calculated from 

the corresponding Mie efficiencies and the size distribution [28] 

00 

f3j = 0.25;r \d^Qj {D)N{D)dD (4.43) 

0 

where D is the diameter of the drop and N(D) is the drop size distribution for the 
rain/fog/haze. 

4.7 RAIN 

Attenuation due to rain is an important problem for infrared communication at 
higher frequencies. In monitoring rainfall and predicting the rain attenuation, the 
raindrop-size distribution plays an important role. 

There are many raindrop size distributions proposed from time to time. We 
analyzed the attenuation due to the rain using the raindrop size distributions, which is 
proposed in the literature. 

4.7,1 Size Distributions for the rain 
Laws and Parsons Distribution 

One X)f the earliest distributions on the size of raindrops was given by Laws and 
Parsons (LP). They proposed the median diameter Dn,: 

= 1 . 238 /?“'“ 

where /? is the precipitation rate in millimeter per hour. 



37 


This distribution has been found to be a reasonable choice for a mean drop-size 
spectrum in continental temperate rainfall, though for rain-rates below 35mm/hr. The LP 
drop size distribution density was fitted to the distribution given by Marshall and Palmer 
(discussed in the next section). 

N{D) = 8.0 X 10" exp(4100i?-®-'D) 
where D is in meter and R is in mm/hr 

Two problems were associated with the above equation are shown below: 

1 . Olsen et al. [29] indicated that it is not consistent with the evaluation of the rain 

rate. 

2. Wolf [30] found that it predicts a maximum value for droplets of zero diameter. 
Wolf [30] proposed a new distribution, which fitted more accurately with the LP data and 
is found to be valid for rain rates of the order of 150 mm/hr. 

The new Laws and Parsons distribution proposed is: 

NiD) = 1.98 X 10’ exp(-5.38/?'" '®"D) (4.44) 

where D is in millimeters and R in mm/hr. 

Marshall and Palmer distribution 

Marshall and Palmer [29] proposed a negative exponential form for the average 
spectrum of drop sizes measured in most rains 

N{D) = No exp(-A D), A= (4.45) 

where No, a and P are constants. 

It has been found to be applicable to widespread rain in continental temperate 
climates, although it has a tendency to overestimate the number of small drops. 

The parameters for the above distribution are: 

No = 0. 08 cni^ and A = 41 Rf '' cm'' 
where R is the rain rate in mm/hr. 
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Thunderstorm distribution 

Joss et al. [29] classified the rain fall into three types, namely, drizzle (J-D), 
widespread (J-W) and thunderstorm (J-T). They proposed a negative exponential 
distribution 

N{D) = Nq exp(-A D), A= aR'^ (4.46) 

with the parameters No = 0.014 cm"^ and A = 30 R'° cm''. 

This negative exponential distribution was fitted to the average drop size measured in 
convective rain. 


Drizzle Distribution (J-D) 

This negative exponential distribution is obtained by fitting the average drop size 
spectrum of veiy light widespread rain, or drizzle, composed mostly of small drops. 
Parameters for this distribution are No = 0.3 cm"^ and A = 57 R'° ^' cm'' 

Widespread Distribution (J-W) 

This negative exponential distribution is obtained by fitting the average drop size 
spectrum of very widespread rain. Parameters for this distribution are No = 0.07 cm"^ and 
A = 4.1 R'®^' cm''. This distribution can be directly obtained by Marshall and Palmer 
distribution as the parameters are identical. Consequently, the values of J-W is equal to 
7/8 of those for the MP distribution. 


Weibull Drop size Distribution 

Jiang et al. [31] indicated that the LP distribution’s behavior is incorrect for very 
small drops. They proposed Weibull distribution which incorporate drizzle, widespread 
rain and shower rain cases. 


c D 


N(D) = N,- - 
by b 


Y-' 


exp( 


f-1 

UJ 


(4.47) 


where D is the diameter in mm. No = 1000 m'^, c - 0.95R^ and b — 0.26R^ ‘’^ mm. 


, 0.14 , 


> 0.44 , 


All the raindrop distribution model talked about till now have been based on the 
temperate regions of the world. The need for the drop size distribution for tropical 
regions was realized and some more models were proposed for the same. 
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Log-normal Distribution 

Ajayi and Olsen [32] pointed out that the MP distribution frequently used for 
scattering applications is inadequate for tropical regions. It overestimates the drop number 
density in the small diameter range at all rain rates, and in the large diameter range at all 
rain rates, and in the large diameter range (above the mode) at high rain rates. 

They proposed a log-normal distribution 


NiD) = 


Nr 

crD^fl^ 



2 




(4.48) 


where p, is the mean of In D, a is the standard deviation, and Nr is the total number of 
drops of all sizes. 

The parameters are: 

Nt= 108 

p = -0.195 + 0.199 In R 
a^ = 0.137 -0.013 In R 

Maitra [33] conducted experiments in Calcutta and gave the following parameters for the 
above distribution 

NT = 54R°•^^^ 
p = -0.538 + 0.0174 In R 
= 0.0685 + 0.0755 In R 


Gamma Distribution 

Gamma distribution was originally proposed by Atlas and Ulbrich [31]. 

N{D) = N, exp(-Ai)), 

where No = 64500 R'°^ m'^ mm"' and A = 7.09 R"®^’ mm"'. 

This distribution was found not to be matching with the tropical region distribution. So, 
Ajayi and Olsen [32] proposed a modified gamma distribution for the tropical regions. 

A(Z)) = 0.0701)' “^ exp(-3 86 D"^), R< 15 mm/hr (4.49) 

N{D) = 0.056 exp(-895 D^ '), R > 15 mm/hr (4.50) 

Maitra [33] mentioned in his paper that only these two parameter model is not adequate 
for describing the drop size distribution in tropical regions and has proposed a three 
parameter gamma model. 
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N(D) = exp(-AD), 

= A,. A"^' D" exp(-AD) / r (« + l) (4.51) 

The parameters are : 

Nt = 564 
n + 1 = 5.61 
A = 9.18R'°^*' 

The lognormal gives lower number density for smaller drops than the gamma distribution. 

4.8 FOG 

Fog presents one of the most important problems to the outdoor optical wireless 
links. In predicting the attenuation due to fog, the drop-size distribution plays an 
important role. 


4.8.1 Size distribution of fog 

Deirmendjian [34] proposed a distribution function of the form 

n(r) = Ar exp(-Z>r’' ) (4.52) 

which vanishes at r = 0, oo and is called as modified gamma distribution. The four 
constants A, a, y and b are positive and real, and a is an integer. They are not 
independent of each other, and are related to quantities in the frequency distribution, 
which can be determined by measurements. 


Upon integration over the entire range of radii, we get 


N = A exp(-ir^ )dr 

0 




(4.53) 


V r J 

which shows that the constant A is essentially given by N, the total number of particles 
per unit volume. Furthermore, differentiation of equation 4.53 with respect to r yields, 


~n(r) = Ar‘'~'{a-ybr’')exp{-br^) (4.54) 

dr 

which has three zeros, two of which are at r = 0 (provided a > 1) and at r = oo, 
respectively. If a = 1 , the derivative at the origin is equal to A. The third zero is found by 
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putting the last factor in equation 4.44 equal to zero, which determines the absolute 
maximum of the function, whence 


with 


n{r^) = A r“ exp 


a 


\ rj 

where rc is the mode radius or size of maximum frequency in the distribution, obtainable 
from the particle counts. Thus the constant b is determined by rc provided a and y are 
fixed or otherwise obtained from the shape of the experimental distribution curve. 


Later, Chu and Hogg [35], modified this distribution, and proposed the following 
modified form of the distribution 
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(4.55) 


where 


A = 


kY. 


N 


a + \ 


and 


b = ^ 
r 


V r ) 

N is the number of drops per unit volume, rc is the mode radius of the drops, and a and y, 
parameters which can be adjusted to fit an observed fog particle size distribution. 


In figure 4.7 four plots of the above equation are given for several combinations 
of parameters. Three of these exhibits a slope on the larger drop side of the maximum (r > 
rc). The other combination of parameters (a = 2, y = 3) corresponds to the asymptotic 
distribution of coagulating particles and has a sharp decreasing rate on both sides of the 
maximum number density. It is of course impossible to have a model for every possible 
fog distribution, but using the different combinations of a and y, as shown in the figure 
4.7, almost all the atmospheric fogs can be described. 
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Using these model distributions, extinction coefficients have been plotted with 
respect the radius at maximum number density rc in figure 4.8 as given in [25]. These 
curves merge at somewhat smaller values of rc because the drops larger than rc contribute 
more to the extinction than the smaller ones. 

The various values of the distributions of the fog and haze for a few a and y 
combinations are tabulated below: 


Distribution Type 

Modal radius (pm) 

a 

a 

Y 

b 

Heavy Fog 

10 

0.027 

3 

1 

0.3 

Moderate Fog 

2 

607.5 

6 

1 

3 

Light Fog 

1 

341 

2 

0.5 

4 . 

Haze L (continental) 

0.07 

5.0e6 

2 

0.5 

15.1 

Tabl 

le 4.2 Distribution model for fog 




43 




lUOtUt OF FMTICLCt W<TM WAXIMUM 

MuMwa ocNtiiT iwicaoMst 


(c) 


- 


S.SOM 

«»t.O 





- ,1 . i M im 



(b) 



OF MttTICl.CS »ITH MAXIMUM 
NUMMtt OCNStTV t MiCttOWS » 


{d} 


Figure 4.8 Extinction cross section for different drop radius [35] 

4.9 SMOG 

Another major challenge faced by the outdoor optical link is the attenuation due to 
smog. Smog is nothing but a mixture of smoke and fog, a colloid system in which the 
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disperse phase consists of a mixture of gas and moisture and the dispersion medium is air. 
Major portion of smog-formers come from burning of petroleum-based fuels such as 
gasoline. Other smog-formers, volatile organic compounds, are found in products such as 
paints and solvents. Major smog occurrences are often linked to heavy motor vehicle 
traffic. In short the smog can be defined as the fog having high pollution content. The 
study of smog is divided into two parts. First the effect of the pollution on the link and the 
study on the overall effect i.e., smog. 

Light extinction due to smoke is made up of the following independent components: 

(a) Particles scatter light in a similar way to gases. Scattering is a function of the 
wavelength of light, particle size and the index of refraction for the particles. Particles 
of 0.5pm in diameter are the most efficient at scattering light. 

(b) Scattering by gases (Rayleigh scattering) which is a constant value for a given 
temperature and pressure. Molecules of gas in the atmosphere divert light from a sight 
path. 

(c) Gases absorb light and transform it into molecular energy. Nitrogen dioxide 
concentration is associated with light absorption by gases because absorption by other 
gas molecules is negligible. But for visible and infrared, there is no molecular 
absorption as shown in figure 4.9. 

(d) Absorption of light by particles occurs when black or colored particles transform light 
into heat. Absorption is linked to the concentration of elemental carbon in the 
atmosphere. 

A number of studies have shown that different chemical species have differing effects 
on visibility degradation. Table 4.3 gives the contribution of the species in the 
degradation of the visibility for visible and near infrared wavelengths. These visibility 
effects are normally split amongst six main groups which are: 

- fine ammonium sulphate 

- fine ammonium nitrate 

- elemental carbon 

- remainder of fine particles (coarse) 
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Species 

% contribution to extinction 

Ammonium Sulphate 

20.2 

Ammonium Nitrate 

17.2 

Carbon 

50 

Particulate Mass (Suspended Particles) 

6.6 

Rayleigh Scattering 

6 


Table 4.3 Visibility apportionment due to pollution 


Following assumptions have been made; 

- all components contributing to scattering are included as independent variables 

- the components are externally mixed 

- scattering efficiency for each component is constant 

The particle scattering coefficient can be defined by the sum of the scattering 
coefficient due to sulphates (bso 4 ), nitrates (bNos) and the coarse suspended particles 
(bcoarse)- Thus the particle scattering coefficient can be expressed as: 

Avp ~ ^N(B ^Ciwse (4.56) 

Because certain aerosols, such as sulphates and nitrates, have an affinity for water, their 
scattering characteristics change as a function of relative humidity (RH). Therefore, 
scattering of the so-called hygroscopic species as a function of relative humidity must be 
considered. The adjustment factor for relative humidity is: 

. f{RH) = 0.34 + 0.59(1 /(I - RH)) + 0.09(1 /(I - RH))' 

Each of the particle scattering coefficients can be related to the mass of the components 
using the following relations: 

= 3[Sulfates)f{RH) 

= 2>{Nitrates]f{RH) 

bc,<ats, - 0.6[Coarse Suspended particles] 

The quantities in brackets are the species concentration expressed in pg/m^ The numeric 
coefficients are the dry specific scattering efficiencies (m''/g). The term f(RH) is the 
relative humidity adjustment factor. The extinction coefficients are in mm . 
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Figure 4.9 Absorption spectra of atmosphere and its constituents[36] 


The absorption due to various gases and molecules is almost zero for the visible 
and the near infrared range. The major absorption that takes place is due to ozone which 
does not affect the terrestrial wireless links. 

So, particle absorption (bap) is primarily due to elemental carbon (soot). 

K, = mEC] 

Again, the quantity in bracket is the mass of elemental carbon in pg/m and 10 is the 
extinction efficiencies. 

The third important parameter that affects the extinction is Rayleigh scattering. As 
the particle size of the gases is very small in comparison to wavelength, Rayleigh 
scattering take place [37] 
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a„,=9.807xlO-»xU” ' W 


T= atmospheric temperature [K] depending on height, 
Pr= pressure [mbar] depending on height, 

X, = laser wavelength [nm]. 


The total extinction for smoke can be expressed as 


Pex! ~ ^S()4 ^Na3 ^Coarse ^an 


(4.57) 


The extinction coefficient of smog is the sum of the extinction coefficient of fog 


and smoke. 


In this chapter, the theory of scattering and radiative transfer is discussed at 
length. The distributions of rain, fog and smog have been discussed and their effect on the 
attenuation on of the received optical signal is studied. In the next chapter, the simulation 
of this theoretical study has been given. 
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CHAPTER 5 

SIMULATION STUDY OF THE OUTDOOR OPTICAL WIRELESS 
LINKS IN RAIN, FOG AND SMOG 

5.1 INTRODUCTION 

Atmospheric channel forms a major part of an outdoor link. As discussed in 
chapter 3, it is one of the major design challenges and affects the link in the form of 
scattering, absorption and turbulence. In chapter 4, we discussed the various means by 
which these parameters come into effect, viz. rain, fog, and smoke. In this chapter, the 
simulation studies on the effect of rain, fog, and smog on the outdoor optical wireless 
links are done and the results are discussed. 

Some assumptions were made while studying the parameters. The particle was 
assumed to be spherical in nature. The medium was considered as homogeneous, i.e, the 
refractive index is constant for the entire link length. For the calculations of the various 
cross sections, Mie-Debye theory has been used as it gives an analysis for scattering and 
absorption by spherical particles. This theory does rigorous treatment of elastic wave 
scattering by spherical particles in a surrounding medium. 

5.2 MIE THEORY 

For the calculation of the Mie efficiency the two important parameters to be 
chosen are the particle size and the refractive index. The refractive index varies with 
variations in the wavelength. The value of the refractive index for various wavelengths 
is taken from reference [38]. Similarly we considered the size of the rain drop to be in 
the range of 0. 1 - 3.6 mm and for fog in the range of 0.0 1 — 40 pm. 

We calculated the extinction efficiency using the equation 4.40 and plotted the 
same as a function of the drop radius in figure 5.1. We see that as the drop radius 
increases, the extinction cross-section approaches the value of 2. Moreover, we find that 
the extinction cross section for higher wavelengths remains substantially smaller than 
the lower wavelengths for radius less than 1 pm. For example, the extinction cross 
section of 3.6 pm is less than that of 0.65 pm only when the radius of water sphere is 
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less than 1 fxm. We see that at 10.0 |xni the first maximum of the extinction cross section 
occurs at a radius of about 25 |xm. 
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Figure 5. 1 Variation of Extinction Cross Section as a function of radius 
for different wavelengths 


The scattering patterns of the electromagnetic waves by spherical particles were 
studied and it was found that as the particle size increases the scattering decreases. The 
scattering patterns for fog and rain were calculated using equation 4.39 and are plotted 
in figure 5.2. From this figure we see that the attenuation due to fog will be much more 
than that of rain. The reason behind this is that for fog reflection and refraction is 
dominant, whereas in the case of rain, diffraction plays a major role. As the particle size 
increases for a fixed wavelength, the diffraction pattern becomes compressed more and 
more into a narrow but very intense lobe around the forward direction, 6 = 0 . Hence for 
such narrow patterns the scattering pattern will be primarily due to diffraction, and for 
broader beams reflection and refraction will play a major role. 
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90 8 



(a) 


90 1000 



(b) 

Figure 5.2 Scattering patterns of electromagnetic wave with spherical particle for 


wavelength of 670 nm (a) Fog (b) Rain 


The scattering and absorption cross sections are calculated using equation 4.4 1 
and 4.42, respectively. The absorption spectrum of the atmosphere for the wavelength 
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range of 0.5—2 |im is studied and is plotted in figure 5.3. It can be seen from the graph 
that for absorption coefficient is highly dependent on the wavelength. For few 
wavelengths, the value of the absorption coefficient is almost zero and hence those 
wavelengths are suitable for the transmission purposes. 



Figure 5.3 Absorption Coefficient as a function of wavelength 
The attenuation for the given absorption coefficient can be calculated from 
equation 4.1 5. This is plotted in figure 5.4. It is clear from the graph that the attenuation 
closely follows the absorption coefficient curve and increases as the absorption 


coefficient increases. 



Figure 5.4 Attenuation for a link distance of 1 km as a function of wavelength 
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The wavelength dependence of the link can be seen from the figure 5.5. We find 
that the absorption for 670 nm is almost zero and that for 10 jim approaches one as the 
radius of the drop increases. As the extinction cross section is the sum of the scattering 
and absorption cross sections, we can conclude that for 670 nm the loss is due to 
scattering only. However, at 10 pm as the radius of the drop increases, the loss is 
contributed equally by absorption and scattering. 




Figure 5.5 Mie Cross sections for wavelength (a) 670 nm (b) 10 pm 
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Figure 5.7 Drop size distribution for different ro- . • 

^''rerent ram rates in mm/hr 

Palmer distribution 


using Marshal and 



Figure 5.8 Drop size distribution for different rain rates in mm/hr 
distribution 


using Drizzle ‘JD’ 


T e „ex, d,se„bu..o„ studied is ,he ‘Drizzle'-JD distribution. E,ua,io„ 4.46 is used to 
caiculate the drop size distribution for diffe^t rain rates and this in plotted in 
fisures.g. This distribution is applicable for vero light widespread rain or drizzle 
contposed mostly of small drops and is valid for teu.perate climates. It can be seen from 
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the graph that the maximum drop size possible for the distribution is only till 1.5 mm. 
The pattern of the distribution resembles the Marshal and Palmer distribution but the 
concentration of the drops is very less for a given volume. 

The next distribution studied is Weibull distribution for which equation 4.47 is used to 
calculate the drop size distribution for different rain rates. We show this result in 
figure 5.9. This distribution is applicable to countries having climate similar to Japan. 
From the graph, we see that as the rain rate increases the concentration for bigger 
particle is more than the particles having smaller diameter. Also, as the rain rate 
increases the diameter of the particle increases. 


X 10'® 



Figure 5.9 Drop size distribution for different rain rates in mm/hr using Weibull 
distribution 

We also studied the gamma distribution and the Log-normal distribution to calculate the 
drop size distribution for different rain rates. Results for the gamma distribution are 
plotted using equation 4.48 in figure 5.10. Log-normal distribution results made use of 
equation 4.49 and are plotted in figure 5.1 1. This distribution is applicable for countries 
having climate similar to India. From the graph, we see that as the rain rate increases the 
concentration of the smaller diameter drop also increases. Also for gamma distribution, 
as the rain rate increases the median diameter of the particle decreases from 0.5 mm to a 
very small value. 
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Figure 5.10 Drop size distribution for different rain rates using Gamma distribution 


X 10‘* 



Figure 5.1 1 Drop size distribution for different rain rates using Log-normal distribution 

After studying all the rain drop distributions, we can generalize the distributions on 
certain points. The contribution of larger rain drops increases with the rainfall rate. 
Regardless of the rainfall rate, most of the contribution of raindrops comes from drops 
of about 0.5 mm diameter. This indicates that the contribution of larger drops is small. 
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5.3,2 Attenuation Characteristics 

Extinction coefficients at different wavelengths as a function of the rain drop 
diameter are plotted in figure 5.12. It is clear from the figure that the scattering due to 
rain is geometric scattering and hence is independent of wavelength. 
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Figure 5.12 Extinction coefficients for different wavelengths as a function of the rain 
drop diameter 


The attenuation coefficient as a function of the rain rate was calculated using 
equations 4.15 and 4.43 and is plotted in figures 5.13. It is found that the attenuation 
increases with the increase in the rain rate. We also see that the attenuation is 
independent of the wavelength employed. 

The attenuation for different rain rates is calculated using equation 4.43 for the 
different drop size distributions. These results are compared in figure 5.14. It can be 
seen from the above plot that for rain rates less than 20 mm/hr almost all the 
distributions show the same attenuation. But as the rain rate increases beyond 20 mm/hr, 
the attenuation calculated for different distributions differs. The lognormal distribution 
is found to overestimate the value of the attenuation at high rain rates and hence is not 
suitable for rain rates greater than 50 mm/hr. Gamma distribution which is valid for 
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Indian sub continent is closely followed by Marshall and Palmer and Weibull 
Distribution. 
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Figure 5.13 Attenuation as a function of rain rate 
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Figure 5.14 Attenuation for different distributions as a function of rain rate 
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5.3.3 Verification of Simulation Result 


We give below the weather parameters from reference [39] for the Huntsville 
Alabama area on July 04 2001. These data will be applied to the theoretical rainfall 
model developed to predict visibility. 

The Laws and Parsons distribution is employed to calculate the attenuation and 
the visibility. 


Time 

Visibility 

Precipitation 

Events 

6:53 PM 

5.0 miles / 8.0 kilometers 

0.02 in / 0.1 cm 

Rain 

8:24 PM 

3.0 miles / 4.8 kilometers 

0.10 in/ 0.3 cm 


8:53 PM 

8.0 miles / 12.9 kilometers 

0.11 in/ 0.3 cm 

Rain 

9:53 PM 

10.0 miles / 16.1 kilometers 

0.06 In / 0.2 cm 

Rain 


Table 5.1 Weather parameters for Huntsville Alabama on July 04, 2001 [39] 


The first simulation was done using the first set of data, i.e., for visibility of 
8.0 km. It is found that the attenuation calculation based on the visibility is greater than 
the one based on drop size distribution. The value of visibility as found from the 
calculation is 10.2 km. 
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Figure 5.15 Received Power as a function of distance 
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The second simulation was done using the second set of data, i.e., for visibility 
of 4.8 km. Here also we find that the attenuation calculation based on the visibility is 
greater than the one based on drop size distribution. The value of visibility as found 
from the model is 5.03 km. 
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Figure 5.16 Received Power as a function of distance 


5.4 FOG 

Fog presents one of the major obstacles for outdoor optical communication. In 
this section, the distribution models of fog and their attenuation characteristics are 
studied. 

5.4.1 Drop Size Distributions 

Drop size distribution for fog changes as the visibility condition due to fog 
changes. The distributions are calculated using Table 4.2. When visibility is less than 
100 m, heavy fog distribution is used. Heavy fog contains particles with radius 0.1 to 40 
pm. 
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Figure 5.17 Particle size distribution for Heavy fog 

When the visibility is less than 500 m, moderate fog distribution is used. Moderate fog 
contains particles with radius 0.01 to 10 pm. 



Figure 5.18 Particle size distribution for Moderate fog 

When the visibility is less than 2 km, light fog distribution is used. Light fog contains 
particles of radius 0.01 to 5 pm. 
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Figure 5,19 Particle size distribution for Light fog 

When the visibility is greater than 2 km, haze distribution is used. Haze contains 
particles with radius 0.01 to 1 pm. 


450 
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Radius (um) 

Figure 5.20 Particle size distribution for Haze 


From the above graphs, it is clear that as the visibility decreases the radius of the fog 
sphere as well as the median radius increases. 
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5.4.2 Attenuation Characteristics 

The extinction coefficients for different wavelengths as a function of the particle 
size are calculated using equation 4.40 and plotted in figure 5.21. We see that as the 
wavelength increases, the extinction efficiency decreases. Also, the peak of the 
extinction efficiency shifts from the lower to higher radii as wavelength increases. 

60 r 



Radius (um) 


Figure 5.21 Extinction coefficients for different wavelengths as a function of the fog 
drop diameter 



Distance (km) 

Figure 5.22 Attenuation as a function of distance for heavy fog 
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Using equation 4.15 we can arrive at the attenuation coefficient for heavy fog at 
different wavelengths. The attenuation for 670nm is found to be 197.9 dB/km and that 
for the 1550nm is 247.46 dB/km. We find that in heavy fogs 1550 nm gets attenuated 
quickly in comparison to 670 nm. The reason for variations in the attenuation in heavy 
fog with wavelength may be the corresponding variations in the extinction coefficient. 


Similarly for moderate fog we can arrive at the attenuation coefficient at different 
wavelengths using equation 4.15. The attenuation for 670nm is found to be 122.5 
dB/km and that for the 1550nm is 125.2 dB/km. The received power for moderate fog 
with distance is shown in figure 5.33. We find that in moderate fogs, the attenuation for 
1550 nm and670 nm is same. 
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Figure 5.23 Attenuation as a function of distance for moderate fog 


A similar approach is used to arrive at the attenuation coefficient for haze. Using 
equation 4.15 we can calculate the attenuation coefficient for haze at different 
wavelengths. The attenuation values for 670nm is found to be 0.18 dB/km and that for 
the 1550nm is 0.06 dB/km The attenuation of the power due to haze with distance is 
plotted in figure 5.24. The haze distribution shows an increase in the attenuation at 670 
nm compared to 1 550 nm. The calculated attenuation for the haze distributions are very 
low and are closer to values expected for clear air. 
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Figure 5.24 Attenuation as a function of distance for haze 


5.5 SMOG 

Smog is another main parameter which can affect the performance of the OWC 
link, especially in places with heavy air pollution. Smog comprises of smoke and fog. In 
this section, the attenuation of the optical power in presence of smog will be discussed. 
The study will be divided into two parts. First the attenuation characteristics in presence 
of smoke will be studied and then this will be combined with the characteristics of fog, 
to get the attenuation due to smog. 

5.5.1 Attenuation Characteristics 

The extinction coefficient for the given concentration of particular matter and 
gas in smoke is calculated using the equation 4.57. Various concentration values for the 
polluted air are taken from the ref [40]. Attenuation over a link distance of 6 km is 
calculated using equation 4.6. 
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Figure 5.25 Attenuation as a function of distance for smoke 

It is seen that even on the days having visibility as high as 200 km in clear air, this 
figure will drop to a few km with smoke (or polluted air). The value of the extinction 
coefficient for smoke is 0.6 km''. Considering dense fog alone the extinction coefficient 
value is as high as 28. Thus the addition of smoke to this condition does not worsen the 
situation too much. However, in light fog and haze where the extinction coefficients are 
typically 4 km'' and 0.05 km'', respectively, addition of smoke leads to very high losses 
over the link distance. The attenuation with distance for different types of fog and smog 
are plotted in figure 5.26 to 5.28. 



Figure 5.26 Attenuation with distance for smog (heavy fog with smoke) and heavy fog 
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Fig 5.27 Attenuation with distance for smog (light fog and smoke) and light fog 
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Figure 5.28 Attenuation with distance for smog(haze and smoke) and haze 


Table 5.2 summarizes the above results and compares the attenuation of power 
over the link for various weather conditions. It can be seen from the table that smog 
affects the link most adversely. If the fog is very dense, then the effect of pollution on 
the link is not very prominent, but for light fog and haze the pollution affects the link 
severely. We also see that the effect of rain is less in comparison to fog and smog. The 
attenuation faced for a rain rate of lOOmm/hr is 1 8 dB/km as compared to an attenuation 
of 206.06 dB/km faced by the link in heavy fog. 
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Link Distance 

Weather 


Attenuation (dB) 

500 m 

Rain 

10 mm/hr 

i.74 



50 mm/hr 

5.53 



1 00 mm/hr 

9.11 


Fog 

Heavy 

103.03 



Moderate 

19.61 



Light 

10.17 


Smog 

Heavy 

104.99 



Moderate 

21.04 



Light 

11.61 

ibOOm 

Rain 

1 0 mm/hr 

3.48 



50 mm/hr 

11.07 



100 mm/hr 

18.22 


Fog 

Heavy 

206.06 



Moderate 

39.21 




20.35 


Smog 

Heavy 

209.98 



Moderate 

42.08 



Light 

2122 

2000 m 

Rain 

10 mm/hr 

6.97 



50 mm/hr 

23.35 



100 mm/hr 

40.71 


Fog 

Heavy 

412.12 



Moderate 

78.44 




40.71 


Smog 


419.97 



Moderate 

84.17 



Light 

46.44 


Table 5.2 Attenuation for different link distances for wavelength = 750 nm and 

transmission power = 1 mW 
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5.6 SIMULATION PACKAGE 

A simulation package has been designed which can perform all the conditions 
listed above. It can calculate the attenuation for a given link in various weather 
conditions. User has to provide some basic information, such as, the transmitter power, 
the wavelength employed, link distance. Given these conditions, the software will 
estimate the output power at the receiving end and will also plot various atmospheric 
effects on the received output power. The GUIs of the simulation package are presented 
below: 
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Figure 5.29 GUI of Simulation package 
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Table 5.3 compares the attenuation over the link under various weather conditions for 
the wavelengths 785 nm and 1550 nm. These results have been generated using the 
simulation package. 


Wavelength 

Weather Conditions 

_ ____ 

Attenuation (dB) 


Rain 

Fog 

Smoke 


785 nm 

Heavy 

- 

- 

18.22 


Drizzle 

- 

- 

1.09 


- 

Heavy 


207.12 


- 

Light 

- 

17.58 


- 

Heavy 


209.98 


- 

Light 


20.02 

1550 nm 

Heavy 


- 

18.22 


Drizzle 

- 

- 

1.09 


- 

Heavy 

- 

247.56 


- 

Light 

1 

26.94 


- 

Heavy 


250.42 


- 

Light 


29.81 

Heavy Rain (Rain rate = 100 mm/hr). Drizzle (Rain rate = 2 mm/1 

hr). 


Heavy Fog (Visibility = 90 m). Light Fog (Visibility = 1 .5 km) 

Table 5.3 Attenuation for a link distance of 1 km and transmission power of 1 m W for 
various weather conditions using the simulation package 
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CHAPTER 6 

CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 


The aim of the present work was to study the effect of atmospheric conditions 
like rain, fog and smog on the performance of outdoor optical wireless 
communication systems. Various weather conditions were studied theoretically and 
their effects on OWC links were simulated. The main results can be summarized as 
below: 


• Smog affects the link most severely. Attenuation due to fog is more than rain. 
In case of fog, reflection and refraction are dominant, whereas in the case of 
rain, diffraction plays a major role. 

• The extinction coefficient for fog varies with the size of the drop but for rain it 
approaches a constant value of 2. For lower wavelengths, the extinction 
coefficient consists of only scattering whereas at higher wavelengths, both 
scattering and absorption plays a dominant role in the calculation of extinction 
coefficient. 

• For all the rain drop size distribution models, the contribution of larger rain 
drops increases with the rainfall rate. Regardless of the rainfall rate, most of 
the contribution of raindrops comes from drops of about 0.5 mm diameter. 
This indicates that the contribution of larger drops is small. The mean of 
raindrop diameters increase with the rainfall rate. 

• The attenuation increases almost linearly with the rain rate and is independent 
of the wavelength employed. The attenuation can be as high as 1 8 dB/km for a 
rain rate of 1 00 mm/hr. 

• In heavy fog conditions, attenuation is very high and hence the operational 
link distance gets reduced to a few meters. 
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• The 1550 nm wavelength light gets attenuated more quickly than the 670 nm 
one in case of heavy fog due to the high value of extinction coefficient. But for 
haze, 670 nm gets attenuated more quickly than 1 550 nm. 

• In very heavy fog, the effect of smoke is not prominent, but for light fog and 
haze conditions, smoke in the air leads to very high losses over the link. 

• While selecting a wavelength the following points should be considered along 
with the weather conditions of the place where the link is to be established. 

a) Attenuation due to rain is independent of wavelength 

b) Attenuation due to heavy fog is high at 1 550 nm. 

c) Attenuation due to light fog and haze is low for 1550 nm. 

6. i SUGGESTION FOR FURTHER WORK 

• This study was basically done for a homogeneous medium. It can be extended 
to incorporate the effects of random non-homogeneous medium, i.e., 
turbulence. 

• In this thesis, only the attenuation of the link is studied. The effect on bit error 
rate and the maximum operable data rate can also be studied. 

• The results of the thesis may be verified by setting up an outdoor optical 
wireless link of length, say 1km, between two buildings and measuring the 
link loss under various weather conditions. 

• The rain, fog, and smog drop-size distribution for the city like Kanpur can be 
studied, which will give an understanding on the various kinds of atmospheric 
losses faced by an optical link. 

• In a country such as India fine dust is a problem during the dry months. The 
effect of dust may also be studied. 
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